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ABSTRACT 

We present an analytic model for blue straggler formation in globular clusters. We 
assume that blue stragglers are formed only through stellar collisions and binary star 
evolution, and co mpare our predictions to obs erved blue straggler numbers taken from 
the catalogue of iLeigh. Sills fc Kniggd ()201lD . We can summarize our key results as 
follows: (1) Binary star evolution consistently dominates blue straggler production in 
all our best-fitting models. (2) In order to account f or the observed sub -linear depen- 
dence of blue straggler numbers on the core masses l|Knigge. Leigh fc Sill s 20091 the 
core binary fraction must be inversely proportional to the total cluster luminosity and 
should always exceed at least a few percent. (3) In at least some clusters, blue strag- 
gler formation must be enhanced by dynamical encounters (either via direct collisions 
or by stimulating mass-transfer to occur by altering the distribution of binary orbital 
parameters) relative to what is expected by assuming a simple population of binaries 
evolving in isolation. (4) The agreement between the predictions of our model and 
the observations can be improved by including blue stragglers that form outside the 
core but later migrate in due to dynamical friction. (5) Longer blue straggler lifetimes 
are preferred in models that include blue stragglers formed outside the core since this 
increases the fraction that will have sufficient time to migrate in via dynamical friction. 
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1 INTRODUCTION 

Commonly found in both open and globular clusters 
(GCs), blue stragglers (BSs) appear as an extension of the 
main-sequence (MS) in cluster colour-magnitude diagrams 
(CMDs), occuyping the region that is just brig hter and 
bluer than the main-sequence turn-off (MSTO) (jSandagd 
1 19531 ). BSs are thought to be produced via the addition 
of hydrogen to low- mass MS stars (e.g. ISills et all l200ll : 
iLombardi et al.ll2002t ). This can occur via multiple channels, 
most of which involve the mergers of low-mass MS stars 
since a significant amount of mass is typically required to 
reproduce the obse rved locations of BSs in CMDs (e.g. 
ISills fc Bailvnlll999l ). Stars in close binaries can merge if 
enough orbital angular momentum is lost, which can be 
mediated by dynamical interactions with other stars, mag- 
netized stellar w inds, tidal dissipation or even an outer triple 
companion fe.g. [Leonard fc Linne ll"l992': 'Li fc Zhang"2006l: 
' I2OO9I : lOervisoglu. Tout fc Ibanoriul 



bOld l. Alternatively, MS stars can collide directly, although 
this is also thoug ht to usually b e mediated by multiple 
star systems 



Fregeau et al.l 



(e.g. [Leonard! '1989: Le onard fc Livid 1 19951 : 



20041 : iLeigh fc Sills 2OI1I). First proposed by 



iPerets fc Fabrvckv 
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McCreal (1964), BSs have also been hypothesized to form 
by mass-transfer from an evolving primary onto a normal 
MS companion via Roche lobe overflow. 

Despite numerous formation mechanisms having been 
proposed, a satisfactory explanation to account for the pres- 
ence of BSs in star clusters eludes us still. Whatever the 
dominant BS formation mechanism(s) operating in dense 
clusters, it is now thought to somehow involve multiple 
star systems. This was shown to be the case in even 
the dense cores of GCs jLeigh. Sills fc Kniggd |2007| . l2008l : 
iKnigge. Leigh fc Silld[200^ where collisions between single 
stars are thought to occur frequently (|Leonardl [19891 '). In 
iKnigge. Leigh fc Silldl|2009l l. we showed that the numbers of 
BSs in the cores of a large sample of GCs correlate with the 
core masses. We argued that our results are consistent with 
what is expected if BSs are descended from binary stars since 
this would imply a dependence of the form Nbs ~ ftMcore, 
where Nbs is the number of BSs in the core, ft is the binary 
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fraction in the core an d Mcors is the total stell ar mass con- 
tained within the core. iMathieu fc Gelleij (|2009l ') also showed 
that at least 76% of the BSs in the old open cluster NGC 188 
have binary companions. Although the nature of these com- 
panions remains unknown, it is clear that binaries played a 
role in the formation of these BSs. 

Blue stragglers are typically concentrated in the dense 
cores of globular clusters where the high stellar densities 
should result in a higher rate of stellar encounters (e.g. 
[Leonard! [T989h . Whether or not this fact is directly related 
to BS formation remains unclear, since mass segregation 
also acts to migrate BSs (or their progenitors) into the core 
(e.g. ISaviane et al.lll998l : Icuhathakurta et al.lll998l '). Addi- 
tionally, numerous BSs have b een observed in more spa rsely 
populated open clusters (e.g. lAndrievskv et all 120001 ) and 
the fields of GCs where collisions are much less likely to oc- 
cur and mass-transfer within binary s ystems is thought to be 
a more likely formation scenario (e.g. iMapelli et alll2004l ). 

Several studies have provided evidence t hat BSs show 
a bimodal s patial distribution i n some GCs (|Ferraro et al.l 
1 19971 . Il999l : iLanzoni et "aIll2007^ . In these clusters, the BS 
numbers are the highest in the central cluster regions and 
decrease with increasing distance from the cluster centre un- 
til a second rise occurs in the cluster outskirts. This drop 
in BS numbers at intermediate cluster radii is often re- 
ferred to as the "zone of avoidance" . Some authors have sug- 
gested that it is the result of two separate formation mech- 
anisms occurring in the inner and outer regions of the clus- 
ter, with mass-transfer in primordial binaries dominating 
in the latter and stellar collisions domin ating in the former 
l|Ferraro et al.ll2004l : iMapelli et all 120061 ). Conversely, mass 
segregation could also give rise to a "zone of avoidance" for 
BSs if the time-scale for dynamical friction exceeds the av- 
erage BS lifetime in only the out skirts of GCs that exhibit 
this radial trend (e.g. iLeigh. Sills & Knigge 2011). 

Dynamical interactions occur frequently enough in 
dense clusters that they are expected to be at leas t partly 
responsible for the obser ved properties of BSs (e.g. IStrvkeil 
ll993l : lLeigh fc Sills|[2oTll ). It follows that the current prop- 
erties of BS populations should reflect the dynamical histo- 
ries of their host clusters. As a result, BSs could provide an 
indirect means of probing the physical processes that drive 
star cluster evolution (e.g . iHeggie fc HutllioO^ : iHurlev et al.l 
l2005l : lLeigh fc Sillsll201ll ). 

In this paper, our goal is to constrain the dominant BS 
formation mechanism(s) operating in the dense cores of GCs 
by analyzing the principal processes thought to influence 
their production. To this end, we use an analytic treatment 
to obtain predictions for the number of BSs expected to be 
found within one core radius of the cluster centre at the 
current cluster age. Predicted numbers for the core are cal- 
culated for a range of free parameters, and then compared to 
the observed numbers in order to find the best-fitting model 
parameters. In this way, we are able to quantify the degree to 
which each of the considered formation mechanisms should 
contribute to the total predicted numbers in order to best 
reproduce the observations. 

In Section[2l we describe the BS catalogue used for com- 
parison to our model predictions. In Section O we present 
our analytic model for BS formation as well as the statisti- 
cal technique we have developed to compare its predictions 
to the observations. These predictions are then compared to 



the observations in Section [4] for a range of model parame- 
ters. In Section [5l we discuss the implications of our results 
for BS formation, as well as the role played by the cluster 
dynamics in shaping the current properties of BS popula- 
tions. 



2 THE DATA 



taken from 
paper, we pre- 



The data used i n this study was 
iLeigh. Sills fc Kniggd (|201ll ). In that 
sented a catalogue for blue straggler, red giant branch 
(RGB), horizontal branch (HB) and main-sequence turn-off 
stars obtained from the colour-magnitude diagrams of 35 
Milky W ay GCs taken from th e ACS Survey for Globular 
Clusters (|Saraiedini et al.ll2007l ). The ACS Survey provides 
unpecedented deep photometry in the F606W (~ V) and 
F814W (~ I) filters that extends reliably from the HB all 
the way down to about 7 magnitudes below the MSTO. 
The clusters in our sample span a range of total masses (by 
nea rly 3 orders of magnitude) and central concentrations 
(Ha rris et al.|[l99^ . We have confirmed that the photom- 
etry is nearly complete in the BS region of the CMD for 
every cluster in our sample. This was done using t he res ults 
of artificial star tests taken from I Anderson et all (|2008l ). 

Each cluster was centred in the ACS field, which ex- 
tends out to several core radii from the cluster centre in 
most of the cluster s in our sample. Only the c ore popu- 
lations provided in iLeigh. Sills fc Knigge! l|201l! ) are used 
in this paper. We have taken estimates for the core radii 
and central luminosity densitie s for the clusters in our 
sample from [Harris et al.l lll996|). whereas cen tral velocity 
dispersions were taken from " IWebbinkI l|l985h . Estimates 
for the total stellar mass contained within the core were 
obtained from single - mass King models, as described in 
iLeigh. Sills fc Kniggd (|201ll ). All of the clusters m our sam- 
ple were chosen to be non-post-core collapse, and have sur- 
face brightness profiles that provide good fits to our King 
models. 



3 METHOD 

In this section, we present our model and outline our as- 
sumptions. We also present the statistical technique used to 
compare the observed number counts to our model predic- 
tions in order to identify the best-fitting model parameters. 



3.1 Model 

Consider a GC core that is home to ^bs,o BSs at some time 
t = to. At a specified time in the future, the number of BSs 
in the core can be approximated by: 



Ne 



NbS,0 + N^oll + A^Mn + Ni„ - Nout ~ iVe. 



(1) 



where NcoH is the number of BSs formed from collisions 
during single-single (1-1-1), single-binary (H-2) and binary- 
binary (2-1-2) encounters, Nj,i„ is the number formed from 
binary evolution (either partial mass-transfer between the 
binary components or their complete coalescence), Ni„ is 
the number of BSs that migrate into the core due to dy- 
namical friction, Nouj is the number that migrate out of the 
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core via kicks experienced during dynamical encounters, and 
Net, is the number of BSs that have evolved away from being 
brighter and bluer than the MSTO in the cluster CMD due 
to stellar evolution. 

We adopt an average stellar mass of m = 0.65Mq and 
an average BS mass of niBS = '2m = 1.3AfQ. The mass of 
a BS can provide a rough guide to its lifetime, although a 
range of lifetimes are still possible for any giv en mass. For 
instance, Sandguist. Bolte fc Hernguistl (|l997l ) showed that 
a 1.3 Mq blue straggler will have a lifetime of around 0.78 
Gyrs in unmixed models, or 1.57 Gyrs i n completel y mixe d 
models. Combined w ith t he results of [ Sills et all l|200ll ). 
iLombardi et"al] (|2002l ) and lGlebbeek fc Polj l|2o'oi ), we ex- 
pect a lifetime in the range 1-5 Gyrs for a 1.3 Mq BS. As 
a first approximation, we choose a lik ely value of tbs = 1.5 
Gyrs for the average BS lifetime (e.g. lSills etmbOOll ). The 
effects had on our results by changing our assumption for 
the average BS lifetime will be explored in Section U and 
discussed in Section (5) 

We consider only the last tbs years. This is because 
we are comparing our model predictions to current observa- 
tions of BS populations, so that we are only concerned with 
those BSs formed within the last few Gyrs. Any BSs formed 
before this would have evolved away from being brighter 
and bluer than the MSTO by the current cluster age. Con- 
sequently, we set N_BS,o = Ne« in Equation [1] We further 
assume that all central cluster parameters have not changed 
in the last tbs years, including the central velocity disper- 
sion, the central luminosity density, the core radius and the 
core binary fraction. It follows that the rate of BS forma- 
tion is constant for the time-scale of interest. This time-scale 
is comparable to the half-mass relaxation time but much 
longer than the central rela xation time for the majority of 
the clusters in our sample ([Harris et al.l Il996l) . This sug- 
gests that core parameters such as the central density and 
the core radius will typically change in a time tbs since the 
time-scale on which these paramet ers vary is the central re- 
laxation time l|Heggie fc Hutll2003l ). Therefore, our assump- 
tion of constant rates and cluster parameters is not strictly 
correct, however it provides a suitable starting point for our 
model. We will discuss the implications of our assumption of 
time-independent cluster properties and rates in Section [S] 

In the following sections, we discuss each of the remain- 
ing terms in Equation [T] 

3.1.1 Stellar Collisions 

We can approximate the number of BSs formed in the last 
Tbs years from collisions during dynamical encounters as: 



Ncoll ~ /l+lA^H-l + /l+2Afl+2 4- f2+2N2+2, 



(2) 



where Ni+i, N1+2 and N2+2 are the number of single- 
single, single-binary and binary-binary encounters, respec- 
tively. The terms fi+i, fi+2 and f2+2 are the fraction of 1-1-1, 
1-1-2 and 2-1-2 encounters, respectively, that will produce a 
BS in the last tbs years. We treat these three variables as 
free parameters since we do not know what fraction of colli- 
sion products will produce BSs (i.e. stars with an appropri- 
ate combination of colour and brightness to end up in the BS 
region of the CMD), nor do we know what fraction of 1-1-2 
and 2-1-2 encounters will result in a stellar collision. Numeri- 
cal scattering experiments have been performed to study the 



outcomes of 1-1-2 and 2-1-2 encounters (e .g. iHut fc Bahcalll 
1 19831 : lMcMillan|[l986l : iFregeau et a"Lll2004 ). however a large 
fraction of the relevant parameter space has yet to be ex- 
plored. 

In terms of the core radius rc (in parsecs), the central 
number density no (in pc~^), the root-mean-square velocity 
Vm (in km s~^), the average stellar mass m (in Mq) and the 
average stellar radius R (in R©), the mean time- scale be- 
tween single-single collisions in the core of a GC is l|Leonardl 
Il989t ): 



n+i = 1.1 X 10"'(l-/f,)-' 



V 5km/s 



0.5Mq\/0.5R, 



(3) 



years 



The additional factor (l-f;,)"^ comes from the fact that we 
are only considering interactions between single stars and 
the central number density of single stars is given by (1- 
f(,)no, where f;, is the binary fraction in the core (i.e. the 
fraction of objects that are binaries). For our chosen mass, 
we assume a correspondin g average s tellar radius using the 
relation M/M© = R/Rq l|lbenlll99ll '). The number of l-fl 
collisions expected to have occurred in the last tbs years is 
then approximated by: 



iVi- 



tbs 

Tl + l ' 



(4) 



The rate of collisions between single stars and binaries, 
as well as between two binary pairs, can be roughly ap- 
proximated in t he same way as for single-sing l e encounters 
jLeonard' 1989; Sigurd sson fc PhinnevI 1 19931 : iBacon et al.l 
[l996.: Fregeau et al. 20 0^1. We adopt the time-scales derived 
m iLeigh fc SillTjioilh for the average times between 1+2 
and 2-1-2 encounters. These are: 

Ti+2 = 3.Axw'f-\i- ur 

'0.5Mq- 



1417 ^ 



lO^pc" 
no 



(5) 



and 



\5km/s/ 



T2+2 = 1.3 X lOVt^' 



years 



i- M 

V 5km/ s / \ m 



)( 



/ lpc\3 flO^ pc 
0.5Mq \ / lAU 



)( 

years. 



r 



(6) 



where a is the average binary semi-major axis in the core in 
AU and we have assumed that the average binary mass is 
equal to twice the average single star mass. The numbers of 
1-1-2 and 2+2 encounters expected to have occurred in the 
last Tbs years are given by, respectively: 

Tbs 



Ni+2 = 



and 



N2 



n+2 

TBS 
T2+2 



(7) 



(8) 



The outcomes of 1-1-2 and 2-|-2 encounters will ulti- 
mately contribute to the evolution of the binary fraction 
in the core. How and with what frequency binary harden- 
ing/softening as well as capture, exchange and ionization 
interactions affect the binary fraction in th e dense cores 
of GCs is currently a subject of debate (e.g. Ilvanova et al.l 
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l2005l : lHurle~ Aarseth fc Sharall2007l ). Observations are also 
lacking for binary fractions in the dense cores of GCs, how- 
ever rough constraints suggest that they range from a few 



to a few tens of a percent (e.g. Rubenstein fc Bailvn 1 19971 : 



ICool et al.|[2003 : ISollima et al.ll2008l : lDavis et aljbooa ). The 



situation is even worse for the distribution of binary orbital 
parameters observed in dense stellar environments. Our best 
constraints come from radial veloci ty surve ys of moderately 
dense open clusters ([Latham. 20051 : iGeller at al. 2009 '). how- 
ever whether or not the properties of the binary populations 
in these clusters should differ significantly from those in the 
much denser cores of GCs is unclear. As an initial assump- 
tion, we assume a time-independent core binary fraction of 
10% for all clusters, and an average semi-major axis of 2 AU. 
This semi-major axis corresponds roughly to the hard-soft 
boundary for most of the clusters in our sample, defined by 
setting the average binary orbital energy equa l to the kinetic 
energ y of an average star in the cluster (e.g. iHeggie fc HutI 
I2OO3I ). We treat both the core binary fraction and the aver- 
age binary semi-major axis as free parameters, and explore 
a range of possibilities using the available observations as a 
guide for realistic values. We will return to these assump- 
tions in Section [5l 



3.1.2 Binary Star Evolution 

Although we do not know the rate of BS formation from 
binary star evolution, we expect a general relation of the 
form Ni,i„ = TBs/Tmt for the number of BSs produced from 
binary mergers in the last tbs years, where Tmt is the av- 
erage time between BS formation events due to binary star 
evolution. We can express the number of BSs formed from 
binary star evolution in the last tbs years as: 



Nbin — fmtfbNcor 



(9) 



where Ncore is the total number of objects (i.e. single and 
binary stars) in the core and fmt is the fraction of binary 
stars that evolved internally to form BSs within the last 
Tbs years. We treat fmt as a free parameter since it is likely 
to depend on the mass-ratio, period and eccentricity distri- 
butions characteristic of the binary populations of evolved 
GC cores, for which data is scarce at best. 



3.1.3 Migration Into and Out of the Core 

Due to the relatively small sizes of the BS populations con- 
sidered, the migration of BSs into or out of the core is an im- 
portant consideration when calculating the predicted num- 
bers. In other words, we are dealing with relatively small 
number statistics and every blue straggler counts. In order 
to approximate the number of stars in the core as a function 
of time, two competing effects need to be taken into account: 
(1) mass stratification/segregation (or, equivalently, dynam- 
ical friction) and (2) kicks experienced during dynamical in- 
teractions. These effects are accounted for with the variables 
Ni„ and Nout in Equation [U respectively. 

Blue st ragglers are among the most massive stars in 
clust e rs (e.g. Shara. Saffer fc Lividll997l : lvan den Berg et al.l 



I2OOII : iMathieu fc Gelleij 120091 ') . so they should typically be 
heavi ly mass segrega,ted relative to oth e r stellar populations 
(e.g. ISpitzeil Il969l : IShara et al.l 1 19951 : iKing. Sosin fc Cooll 



Il995h . The time-scale for this process to occur can 
be approximated using t he dynamical friction time-scale 
(|Binnev fc Tremaind[l987l ): 



a{rf 



4ZnAG2(27r)i/2 mssp(r-)' 



(10) 



where (j{r) and p{r) are, respectively, the velocity dispersion 
and stellar mass density at the given distance from the clus- 
ter centre r. Both (j(r) and p(r) are found from single-mass 
King models l|Sigurdsson fc Phinnevl[l99^ . which are fit to 
ea ch cluster using the c oncentration parameters provided 
in McLa ughlin fc van de n Marol (2005). The Coulomb log- 
arithm is denoted by A, and we adop t a va l ue of In A ~ 10 
throughout this paper (e.g. ISpitzed Il987l : iHeggie fc HutI 
l2003h . If Tdyn > TBS at a given distance from the cluster 
centre, then any BSs formed at this radius in the last tbs 
years will not have had sufficient time to migrate into the 
core by the current cluster age. The maximum radius Tmax 
at which BSs can have formed in the last tbs years and still 
have time to migrate into the core via dynamical friction is 
given by setting Tdyn = tbs - Therefore, Ni„ depends only on 
the number of BSs formed in the last tbs years at a distance 
from the cluster centre smaller than Tmax. 

In order to estimate the contribution to Nss in Equa- 
tion [T] from BSs formed outside the core, we calculate the 
number of BSs formed in radial shells between the cluster 
centre and Tmax. Each shell is taken to be one core radius 
thick, and we calculate the contribution from each formation 
mechanism in every shell. This is done by assuming a con- 
stant (average) density and velocity dispersion in each shell. 
Specifically, we estimated the density and velocity dispersion 
at the half-way point in each shell using our single-mass King 
models, and used these to set average values. The number 
of BSs expected to have migrated into the core within the 
last Tbs years can be written: 

JV 

Nin = ^ (^/l+lAi'(i + i),i -I- fl+2N^i+2),i + /2-l-2iV(2+2),i 

(11) 

where i = 1 refers to the core, i — 2 refers to the shell 
immediately outside the core, etc. and N is taken to be the 
integer nearest to Vmax/rc. We let the terms with N(i+i) 
N{i+2),i, N(2+2),i and N(^tin),i represent the number of BSs 
formed in shell i from single-single collisions, single-binary 
collisions, binary-binary collisions and binary star evolution, 
respectively. The time-scale for dynamical friction in shell 
i is denoted by T(^dyn),ii and the factor (1 - T^dyn), z/tbs) 
is included to account for the fact that we are assuming a 
constant formation rate for BSs, so that not every BS formed 
in shells outside the core will have sufficient time to fall in 
by the current cluster age. 

It is typically the least massive stars that are ejected 
from 1-1-2 and 2-|-2 inter actions as single stars (e.g. 
ISigurdsson fc Phinnevl[l993 ). Combined with conservation 
of momentum, this suggests that BSs are the least likely 
to be ejected from dynamical encounters with velocities 
higher than the central velocity dispersion due to their large 
masses. This has been confirmed by several studies of nu- 
merical scattering experiments (e.g. iHut fc Bahcalll Il983l : 
iFregeau et al.l [2004 ). Based on this, we expect that Nout 
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should be very small and so, as a first approximation, we 
take Nout = 0. However, we also explore the effects of a non- 
zero Nout by assigning a kick velocity to all BSs at birth. If 
dynamical interactions play a role in BS formation, we might 
naturally expect BSs to be imparted a recoil velocity at birth 
(or shortly before) due to momentum conservation. We will 
return to this assumption in Section |4l 



4 RESULTS 

In this section, we present the results of comparing our 
model predictions to the observations. After presenting the 
results for a constant core binary fraction for all clusters, 
we explore the implications of adopting a core binary frac- 
tion that dep e nds o n the cluster luminosity, as reported in 
ISollima et atl (|2007|) andlMilone et all (|2008|). 



3.2 Statistical Comparison with Observations 

Our model contains 4 free parameters, which corre- 
spond to the fraction of outcomes that produce a blue 
straggler for each formation mechanism (H-l collisions, 
1-1-2 collisions, 2-1-2 collisions, and binary star evolution). 
These are the / values described in the previous section: 
f i+i, f 1+2, f 2+2, fmt- We assume that these values are con- 
stant throughout each cluster, and are also constant between 
clusters. By fitting the predictions of our model to the obser- 
vational data, we can determine best-fit values for each of 
these / parameters, and therefore make predictions about 
which blue straggler formation processes are more impor- 
tant. 

In order to determine the best values for these / pa- 
rameters, we need an ap propriate statistical test. For this , 
we follow the approach of lVerbunt. Poolev fc Bassal l|2008h . 
The number of BSs seen in the core of a globular cluster can 
be described by Poisson statistics. In particular, the proba- 
bility of observing A'' sources when fi are expected is: 

P(iV,/i) = ^e-'* (12) 

We can calculate a probability for each cluster, and then cal- 
culate an overall probability P' for the model by multiplying 
the individual P values. We can then vary the / values to 
maximize this value. 

In practice, these P values are typically of order ten per- 
cent per cluster, and with 35 clusters, the value of P' quickly 
becomes extremely small. Therefore we chose to work with a 
modified version of this value: the deviance of our model to 
the saturated model. A saturated model is one in which the 
observed number of sources is exactly equal to the expected 
number in each cluster. In other words, this is the best that 
we can possibly do. However, because of the nature of Pois- 
son statistics, the probability P of such a model (calculated 
by setting N — fi in Equation 1 12|) is not equal to 1, but in 
fact has some smaller value. For the numbers of blue strag- 
glers in our clusters, the P values for the saturated model 
run from 0.044 to 0.149, and the value of P' is 2.08 x lO^^V 

The deviance of any model from the saturated model is 
given by 

D = 2.0(ln(Pl,„„,,^) - ln(P4„d,,)) (13) 

The model which minimizes this quantity will be our best-fit 
model. Ideally, the scaled deviance [D / (Ndata- Nparameters) 
should be equal to 1 for a best fit. Given the simplicity of 
our model, we expect that our values will not provide this 
kind of agreement, and we simply look for the model which 
provides the minimum of the scaled deviance. 



4.1 Initial Assumptions 

The predictions of our model for our initial choice of as- 
sumptions are shown in Figure [1] These numbers are plot- 
ted against the total stellar mass in the core along with the 
number of BSs observed in the core (filled triangles). We 
plot both the total number of BSs predicted to have formed 
within rmax in the last tbs years (N_bs in Equation[T] open 
circles) , as well as the total number formed only in the core 
(Nco!! + Ntin; small filled circles). Upon comparing Nbs 
to the observed number of BSs in the core, the best-fitting 
model parameters predict that most BSs are formed from 
binary star evolution, with a small contribution from 2-1-2 
collisions being needed in order to obtain the best possible 
match to the observations. The ideal contribution from 2-1-2 
collisions constitutes at most a few percent of the predicted 
total for most of the clusters in our sample. Even for our 
best-fitting model parameters, our initial choice of assump- 
tions predicts too few BSs in clusters with small core masses. 

4.2 Binary Fraction 

We tried changing our assumption of a constant f;, for all 
clusters to one for which the core binary fraction varies with 
the total cluster magnitude. First, we adopted a dependence 
of the form: 

/b = O.lSA/v- + 1.07, (14) 

where Mv is the total cluster V magnitude. This relation 
comes from fi tting a line of best-fit to the observations of 
ISollima et al.l |2007) , who studied the binary fractions in a 
sample of 13 low-density GCs (we calculated an average of 
columns 3 and 4 in their Table 3 and used these binary frac- 
tions to obtain Equation I14p . In order to prevent negative 
binary fractions, we impose a minimum binary fraction of 
f^™ = 0.01. In other words, we set ft = if Equation [Til 
gives a binary fraction less than f^™. As before, we adopt 
an average semi-major axis of 2 AU. The results of this com- 
parison are presented in Figure [5] As in Figure [T] both the 
numbers of observed (filled triangles) and predicted (open 
circles) BSs in the core are plotted versus the total stellar 
mass in the core. Once again, the predicted numbers include 
all BSs formed within Xmax in the last tbs years. The best- 
fitting model parameters for this comparison suggest that 
both single-single collisions and binary star evolution are sig- 
nificant contributors to BS formation. Single-single collisions 
contribute up to several tens of a percent of the predicted to- 
tal in several clusters. We obtain a deviance in this case that 
is significantly larger than that obtained for the best-fitting 
model parameters assuming a constant f(, of 10% for all clus- 
ters. Equation 1141 gives higher binary fractions in low-mass 
clusters relative to our initial assumption of a constant f;,. 
This increases the number of BSs formed from binary star 
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Figure 1. The predicted number of BSs plotted versus the to- 
tal stellar mass in the core for the best-fitting model parameters 
found for our initial choice of assumptions. The filled triangles 
correspond to the observed numbers, the open circles to the num- 
ber of BSs predicted to have formed within Vmax, and the small 
filled circles to the number predicted to have formed in only the 
core. The best-fitting model parameters used to calculate the pre- 
dicted numbers are fi+i = 0, fi-i-2 = 0, f2+2 = 7.3 X 10~^, and 
fmt = 1.7 X 10~^. Estimates for the total stellar masses in the 
core were obtained from sing le-mass King models, as described in 
iLeigh. Sills Sz Kniggd ll201lf ). Error bars have been indicated for 
the observed numbers using Poisson statistics. 



evolution and improves the agreement between our model 
predictions and the observation s in low-mass c luster s. This 
is consistent with the results of ISollirna et al.1 (|2008l '). How- 
ever, adopting Equation [14] for f;, also causes our model to 
under-predict the number of BSs in several high-mass clus- 
ters. 

The best fit to the observations is found by adopting the 
relation for fj, provided in Equation[T4]and setting fJJ*'" = 0.1 
(however we note that a comparably good agreement is 
found with a slightly lower fj^"" = 0.05). This improves the 
agreement between our model predictions and the observa- 
tions by increasing the number of BSs formed from binary 
star evolution in massive clusters. The result is a good agree- 
ment between our model predictions and the observations in 
both low- and high- mass cores, as shown in Figure O In this 
case, the best-fitting model parameters yield the lowest de- 
viance of any of the assumptions so far considered. These 
best-fitting values suggest that most BSs are formed from 
binary star evolution, with a small contribution from 2+2 
collisions being needed in order to obtain the best possi- 
ble match to the observations. Similarly to what was found 
for our initial assumptions, the ideal contribution from 2+2 
collisions constitutes at most a few percent of the predicted 
total for most of the clusters in our sample. On the other 
hand, if we change our imposed minimum binary fraction 
to f^™ — 0.05 we find that a non- negligible (i.e. up to a 
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Figure 2. The number of BSs predict e d in t he cluster core using 
the binary fractions of ISoUima et al.l ll2007t) with f™*" = 0.01 
plotted against the total stellar mass contained within the core. 
The symbols used to indicate the observed and predicted numbers 
are the same as in Figure[T] The predicted numbers correspond to 
the best-fitting model parameters, which are fi-|-i = 0.41, fi+2 = 
0, f2-|-2 = 0, and fmt = 1.5 X lO'^. 

few tens of a percent) contribution from single-single colli- 
sions is needed in several clusters to obtain the best possible 
agreement with the observations (which is very nearly as 
good as was found using f^™ = 0.1). All of this shows that, 
although binary star evolution consistently dominates BS 
formation in our best-fitting models, at least some contri- 
bution from collisions (whether it be 1-f 1 or 2+2 collisions, 
or some combination of 1-1-1, 1-1-2 and 2+2 collisions) also 
consistently improves the agreement with the observations. 
Moreover, it is interesting to note that an improved agree- 
ment with the observations could alternatively be obtained 
if we keep 1^*" = 0.01 but all or some of fi+i, fi+2 and 
f2+2 increase with increasing cluster mass. This would also 
serve to improve the agreement at the high-mass end. We 
will return to this in Section [5] 

Finally, we also tri ed adopting the obse rved dependence 
of fj, on yiv reported in lMilone et al.l l|2008l ). who also found 
evidence for an anti-correlation between the core binary frac- 
tion and the total cluster mass. Despite this change, we con- 
sistently find that our results are the same as found when us- 
ing the empirical binary fraction relation provided in Equa- 
tion [H 

4.3 Average BS Lifetime 

We also tried changing our assumption for the average BS 
lifetime. We explored a range of plausible lifetimes based on 
values found throughout the literature. Specifically, we ex- 
plored the range 0.5-5 Gyrs. We find that at the low end of 
this range, our model fits become increasingly poor. This is 
because lower values for tbs correspond to smaller values for 
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Figure 3. The numbe r of BSs predi c ted in the cluster core using 
the binary fractions of lSoIlima et al.l ||2007|) with f^J^'" = 0.1 plot- 
ted against the total stellar mass contained within the core. The 
symbols used to indicate the observed and predicted numbers are 
the same as in Figure [T] The predicted numbers correspond to 
the best-fitting model parameters, which are fi-)-i = 0, fi+2 = 0, 
h+2 = 3.6 X 10-3, and f^t = 1.4 X 10"^. 

i^max and decrease the term (1 - t(^ayn),i/'''Bs) in EquationfTT] 
This reduces the contribution to the total predicted numbers 
from BSs formed outside the core that faU in via dynamical 
friction. Conversely, our model fits improve for tbs > 1-5 
Gyrs since this corresponds to a larger contribution to Nss 
from Ni„. It is important to note, however, that this same 
effect can be had by increasing the number of BSs formed 
outside the core, since this would also serve to increase Ni„ 
in Equation [T] This can be accomplished by, for instance, 
increasing the binary fraction outside the core (which would 
increase the number of BSs formed from binary star evo- 
lution outside the core that migrate in due to dynamical 
friction) relative to inside the core. This seems unlikely, how- 
ever, given that observations of low-density globular clusters 
and open clusters suggest that their bi nary fractions tend t o 
drop off rapidly outside the core (e.g. ISollima et a^l2007^ . 
We will return to this issue in Section [S] 

The best possible match to the observations is found 
by adopting an average BS lifetime of 5 Gyrs along with the 
relation for f^ provided in Equation [T?] with fj^'" = 0.1. The 
predictions of our model are shown in Figure[4]for these best- 
fitting model parameters. As shown, the agreement between 
our model predictions and the observed numbers is excellent. 
The best agreement is found by adopting an average BS 
lifetime of 5 Gyrs, however the agreement is comparably 
excellent down to slightly less than tbs ~ 3 Gyrs. Although 
increasing tbs does contribute to improving the agreement 
between our model predictions and the observations, the 
effect is minor compared to the improvement that can be 
found by changing our assumption for the binary fraction. 
This is apparent upon comparing Figure |4] to Figure |3j 
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Figure 4. The predicted number of BSs plotted versus the to- 
tal stellar mass in the core for the best-fitting model parame- 
ters found using an average BS lifetime of 5 Gyrs and the re- 
lation for the cluster binary fraction provided in Equation 1141 
with fj^'" =0.1. The symbols used to indicate the observed and 
predicted numbers are the same as in Figure [T] The best- fitting 
model parameters used to calculate the predicted numbers are 
fi+i = 0, fi+2 = 0, f2+2 = 1-6 X 10-3, and f™t = 9.9 x 10-^. The 
agreement with the observations is excellent for these best-fitting 
values. 



4.4 Migration 

In order to explore the sensitivity of our results to our as- 
sumption for Tmax , we also tried setting Ni„ equal to the to- 
tal number of BSs expected to form within 10 Vc for all clus- 
ters. For comparison, for an average BS lifetime of tbs = 1.5 
Gyrs, Tmax ranges from 2 - 15 rc for the clusters in our sam- 
ple. Despite implementing this change, our results remained 
the same. This is because the largest contribution to the to- 
tal number of BSs comes from those BSs formed in the first 
few shells immediately outside the core that migrate in due 
to dynamical friction. 

Several GCs have been reported to show evidence for 
a decrease in their binary fr actions with increasing distance 
from the cluster centre (e.g. ISoUima et al.|[2007l : iDavis et al.l 
l2008h . This effect is often significant, with binary fractions 
decreasing by up to a factor of a few within only a few core 
radii from the cluster centre. Based on this, our assump- 
tion that f(, is independent of the distance from the cluster 
centre likely results in an over-estimate of the true binary 
fraction at large cluster radii. In order to quantify the pos- 
sible implications of this for our results, we tried setting fb 
= for all shells outside the core. Although this assumption 
is certainly an under-estimate for the true binary fraction 
outside the core, our results remain the same (albeit the 
agreement with the observations is considerably worse than 
for most of our previous model assumptions). Once again, 
the best-fitting model parameters suggest that most BSs are 
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formed from binary star evolution, with a non-negligible (i.e. 
up to a few tens of a percent in some clusters) contribution 
from binary-binary collisions. Our results indicate that, if 
the binary fraction is negligible outside the core, then the 
contribution from BSs that migrate into the core due to dy- 
namical friction is also negligible. This is because the time 
between 1-1-1 collisions increases rapidly outside the core, 
and every other BS formation mechanism requires binary 
stars to operate. 

We also explored the effects of assuming a non-zero 
value for No„t in Equation [T] by imparting a constant kick 
velocity to all BSs at birth. Using conservation of energy, we 
calculated the cluster radius to which BSs should be kicked 
upon formation, and used the time-scale for dynamical fric- 
tion at the kick radius to calculate the fraction of BSs ex- 
pected to migrate back into the core within a time tbs- 
Regardless of our assumption for the kick velocity, this did 
not improve the deviance for any of our best-fitting model 
parameters. 



4.5 Average Binary Semi-Major Axis 

We investigated the dependence of our results on our as- 
sumption for the average binary semi-major axis. However, 
this had a very small effect on our results. This is because 
only Ni+2 and N2+2 depend on the average semi-major axis, 
and neither of these terms dominated BS production regard- 
less of our model assumptions. Only f2+2 is non-zero for our 
best-fitting models however, as before, it consistently sug- 
gests that far fewer BSs should be formed from 24-2 collisions 
than from binary star evolution. 



5 SUMMARY & DISCUSSION 

In this paper, we have presented an analytic model to inves- 
tigate BS formation in globular clusters. Our model predicts 
the number of BSs in the cluster core at the current cluster 
age by estimating the number that should have either formed 
there from stellar collisions and binary star evolution, or 
migrated in via dynamical friction after forming outside the 
core. We have compared the results of our model predictions 
for a variety of input parameters to observed BS numbers 
in 35 GCs taken from the catalogue of lLeigh. Sills fc Kniggel 
l|201lf ). 

What has our model told us about BS formation in 
dense cluster environments? The agreement between the 
predictions of our model and the observations is excellent 
if we assume that: 

• Binary star evolution dominates BS formation, how- 
ever at least some contribution from 2-1-2 collisions (most of 
which occur in the core) must also be included in the to- 
tal predicted numbers. Although it is clear that including a 
contribution from dynamical encounters gives the best pos- 
sible match to the observations, it is not clear how exactly 
this is accomplished in real star clusters. Does the cluster 
dynamics increase BS numbers via direct collisions? Or do 
dynamical interactions somehow modify primordial binaries 
to initiate more mass-transfer events? We will return to this 
point below. 



• The binary fraction in the core is inversely correlated 
with the total clus ter luminosity , similar to the empirical 
relations found by ISollima et al. I l|2007l ) and iMilone et al.l 
(2008). We also require a minimum core binary fraction of 
5 — 10%. The inverse dependence of ft, on the total clus- 
ter mass contributes to a better agreement with the obser- 
vations at the low-mass end of the distribution of cluster 
masses, whereas the imposed condition that f[^'" = 0.05—0.1 
contributes to improving the agreement at the high-mass 
end. 

• BSs formed outside the core that migrate in by the cur- 
rent cluster age contribute to the total predicted numbers. 

• The average BS lifetime is roughly a few (~ 3-5) Gyrs, 
since this increases the fraction of BSs formed outside the 
core that will have sufficient time to migrate in due to dy- 
namical friction. 



Our model can only provide a reasonable fit to the ob- 
servations for all cluster masses if we assume that the cluster 
binary fraction is inversely proportional to the total cluster 
mass. It is interesting to consider the possibility that such 
an inverse proportionality could arise as a result of the fact 
that the rate of two-body relaxation is also inversely propor- 
tional to the cluster mass (|Spitzerl 19871 ) . Consequently, since 
binaries tend to be the most massive objects in GCs, they 
should quickly migrate into the core in low-mass clusters, 
contri buting to an increase in the cor e binary fraction over 
time (jFregeau. Ivanova fc Rastoll2009l ). This process should 
operate on a considerably longer time-scale in very mas- 
sive GCs since the time-scale for two-body relaxation is very 
long. Mass segregation could then contribute to the observed 
sub-linear dependence of BS numbers on the core masses by 
acting to preferentially migrate the binary star progenitors 
of BSs into the cores of low-mass clusters. This is one ex- 
ample of how a direct link could arise between the observed 
properties of BS populations and the dynamical histories 
of their host clusters. Although this scenario is interesting 
to consider, we cannot rule out the possiblity that an anti- 
correlation between the core binary fraction and the total 
cluster mass could be a primordial property characteristic 
of GCs at birth. 

When interpreting our results, it is important to bear in 
mind that binary star evolution and dynamical interactions 
involving binaries may not always contribute to BS forma- 
tion independently. For example, dynamics could play an 
important role in changing the distribution of binary orbital 
parameters so that mass-transfer occurs more commonly in 
some clusters. One way to perhaps compensate for this effect 
would be to include a factor of 1/a (where a is the average 
binary semi-major axis) in Equation |5] This would serve to 
account for the fact that we might naively expect clusters 
populated by more close binaries to be more likely to have 
a larger fraction of their binary populations undergo mass- 
transfer. This does not, however, guarantee that more BSs 
will form since our poor understanding of binary star evo- 
lution prevents us from being able to predict the outcomes 
of these mass-transfer events, and whether or not they will 
form BSs. Moreover, little is known about the distribution of 
orbital parameters characteristic of the binary populations 
in globular clusters, and how they are typically modified by 
the cluster dynamics. For these reasons, the interpretation 
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of our results must be done with care in order to ensure that 
reliable conclusions can be drawn. 

In general, our results suggest that binary stars play a 
crucial role in BS formation in dense GCs. In order to ob- 
tain the best possible agreement with the observations, an 
enhancement in BS formation from dynamical encounters 
is required in at least some clusters relative to what is ex- 
pected by assuming a simple population of binaries evolving 
in isolation. It is not clear from our results, however, how 
exactly this occurs in real star clusters. Dynamics could en- 
hance BS formation directly by causing stellar collisions, 
or this could also occur indirectly if the cluster dynamics 
somehow induces episodes of mass-transfer by reducing the 
orbital separations of binaries. But in which clusters is BS 
formation the most strongly influenced by the cluster dy- 
namics? Unfortunately, no clear trends have emerged from 
our analysis that provide a straight-forward answer to this 
question. However, our results are consistent with dynamical 
interactions playing a more significant role in more massive 
clusters (although this does not imply that the cluster dy- 
namics does not also contribute in low-mass clusters). This 
could be due to the fact that more ma ssive clusters also tend 
to ha ve higher central densities (e.g. iDiorgovski fc MevlanI 
1 19941 ). and therefore higher collision rates. This picture is, 
broadly speaking, roughly consis tent with the results of 
iDavies. Piotto fc De Angelil (|2004l '). These authors consid- 
ered the observed dependence (or lack thereof) of BS num- 
bers o n the total cluster masses presented in iPiotto et al.l 
and suggested that primordial binary evolution and 
stellar collisions dominate BS production in low- and high- 
mass clusters, respectively. 

Our model neglects the dynamical evolution of GCs and 
the resulting changes to their global properties, including 
the central density, velocity dispersion, core radius and bi- 
nary fraction. As a young cluster evolves, dynamical pro- 
cesses like mass segregation and stellar evaporation tend 
to result in a smaller, denser core. Within a matter of a 
few half-mass rela:xation times, a gravothermal instability 
has set in and the collapse ensues on a time-scal e deter- 
mine d by the rate of heat flow out of the core (e.g. ISpitzej 
1 19871 ). We are focussing on the last tbs years of cluster evo- 
lution, a sufficiently late period in the lives of most GCs that 
gravothermal collapse will have long since taken over as the 
primary driving force affecting the stellar concentration in 
the core. Most of the GCs in ou r sample should currently 
be in a phase of core contraction |Fregeau. Ivanova fc Rasid 
I2OO9I : ICieles. Heggie fc Zhaoll201ll ). and their central den- 
sities and core radii should have been steadily decreasing 
over the last tbs years. Therefore, by using the currently 
observed central cluster parameters and assuming that they 
remained constant over the last tbs years, we have effec- 
tively calculated upper limits for the encounter rates. This 
could suggest that we have over-estimated the importance 
of dynamical interactions for BS formation. On the other 
hand, some theoretical models of GC evolution suggest that 
the hard binary fraction in the core of a dense stellar sys- 
tem will generally i ncrease with time (e.g. lHurlev et al" ]|2005l : 
iFregeau. Ivanovafc Rasio 2009) . This can be understood as 
an imbalance between the migration of binaries into the core 
via mass segregation and the destruction of binaries in the 
core via both dynamical encounters and their internal evolu- 
tion. This could suggest that our estimate for the number of 



BSs formed from binary star evolution should also be taken 
as an upper limit. The key point is that GC evolution can 
act to increase the number of BSs in the cluster core via sev- 
eral different channels. The effects we have discussed should 
typically be sm all, however, since tbs is much shorter than 
the cluster age (|De Angeh et al.ll2005l ) and any changes to 
global cluster properties that occur during this time will of- 
ten be small. 

Our model adopts the same values for all free param- 
eters in all clusters. In particular, this is the case for sev- 
eral global cluster properties, including the average stellar 
mass, the average BS mass, the average BS lifetime, and 
the average binary semi-major axis. With the exception 
of the average stellar mass, there is no conclusive obser- 
vational or theoretical evidence to indicate that these pa- 
rameters should differ from cluster-to-cluster, although we 
cannot rule out this possibility. For instance, the distribu- 
tion of binary orbital parameters could depend on cluster 
properties like the total mass, density or velocity disper- 
sion (e.g. ISigurdsson fc PhinnevI ri993l ). In particular, the 
central velocity dispersion should be highe r in more mas- 
sive GCs (e.g. iDiorgovski fc MevlanI 1 19941 ). which should 
correspond to a smaller binary orbital separation for the 
hard-soft boundary. This could contribute to massive GCs 
tending to have smaller average binary orbital separations 
since soft binaries should not sur vive for long in the dense 
cores of GCs (e.g. iHeggie fc Hutl l2003). In turn, this could 
affect the occurrence of mass-transfer events, or of merg- 
ers during l-|-2 and 2-|-2 encounters. This last point follows 
from the fact that numerical scattering experiments have 
shown that the probability of mergers occurring during 1+2 
and 2+2 interaction s increases with dec reasing binary or- 
bital separation (e.g. iFregeau et al]|2004h . Both the average 
stellar mass and the average BS mass (and hence lifetime) 
co uld also depend on th e tota. 1 clust er mass, as discussed 



■Leigh. Sills fc Kniggd (|2009l ) and iLeigh. Sills fc Kniggd 
(|201lf ). 

We have also neglected to consider the importance 
of triples for BS f ormation throughout our analysis (e.g. 
iPerets fc Fabrvckvl [2009) since we are unaware of any ob- 
servations of triples in GCs in the literature. Interestingly, 
however, our results for binary star evolution can be gener- 
alized to include the internal evolution of triples since they 
should have the same functional dependence on the core 
mass (i.e. Nte oc ftMcore, where Nte is the number of BSs 
formed from triple star evolution and ft is the fraction of 
objects that are triples). 

Finally, our model assumes that several parameters re- 
main constant as a function of the distance from the cluster 
centre, including the binary fraction and the average semi- 
major axis. However, observations of GCs suggest that their 
binary fr actions could fall off rapidly outside the core (e.g. 
Sollima et al. 2007; Davis ct al. 2008). Our results suggest 
that, if the binary fraction is negligible outside the core, then 
the contribution from BSs that migrate into the core due to 
dynamical friction is also negligible. This is because the time 
between 1+1 collisions increases rapidly outside the core, 
and every other BS formation mechanism requires binary 
stars to operate. On the other hand, the presently observed 
binary fraction outside the core could be low as a result 
of binaries having previ ously migrated into the core du e to 
dynamical friction (e.g. IFregeau. Ivanova fc Rasid [20091 ) . If 
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these are the binary star progenitors of the BSs currently 
populating the core, then dynamical friction remains an im- 
portant effect in determining the number of BSs currently 
populating the core. 

Despite all of these simplifying assumptions, wc have 
shown that our model can reproduce the observations with 
remarkable accuracy. Notwithstanding, the effects we have 
discussed could be contributing to cluster-to-cluster differ- 
ences in the observed BS numbers. Our model provides a 
well-suited resource for addressing the role played by these 
effects, however future observations will be needed in order 
to obtain the desired constraints (e.g. binary fractions, dis- 
tributions of binary orbital parameters, etc.). 
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